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Q . ABSTRACT 

^■^ \ Integrated star formation rate (SFR) and specific star formation rate (SFR/M*), derived 

r"| . from the spectroscopic data obtained by SDSS DR4 are used to show that the star formation 

Q_[ activity in galaxies (M r ^ —20.5) found on the outskirts (l-2r 2 oo) of some nearby clusters 

(0.02 ^ z ^ 0.15) is enhanced. By comparing the mean SFR of galaxies in a sample of clus- 
ters with at least one starburst galaxy (log SFR/A/* ^ -10 yr _1 and SFR ^ 10 Moyr^ 1 ) to a 
£^ 1 sample of clusters without such galaxies ('comparison' clusters), we find that despite the ex- 

■ pected decline in the mean SFR of galaxies toward the cluster core, the SFR profile of the two 

samples is different. Compared to the clusters with at least one starburst galaxy on their out- 
skirts, the galaxies in the 'comparison' clusters show a low mean SFR at all radius (^ 3r 2 oo) 
from the cluster centre. Such an increase in the SFR of galaxies is more likely to be seen in 
dynamically unrelaxed (a v > 500 km s _1 ) clusters. It is also evident that these unrelaxed 
clusters are currently being assembled via galaxies falling in through straight filaments, re- 

■ suiting in high velocity dispersions. On the other hand, 'comparison' clusters are more likely 
C^> ' to be fed by relatively low density filaments. We find that the starburst galaxies on the pe- 
riphery of clusters are in an environment of higher local density than other cluster galaxies 

£^ ■ at similar radial distances from the cluster centre. We conclude that a relatively high galaxy 

| density in the infalling regions of clusters promotes interactions amongst galaxies, leading to 

momentary bursts of star formation. Such interactions play a crucial role in exhausting the 
fuel for star formation in a galaxy, before it is expelled due to the environmental processes 
that are operational in the dense interiors of the cluster. 

. Key words: Galaxies: clusters: general; Galaxies: evolution; Cosmology: observations; Cos- 

■ mology: large-scale structure of the Universe; Galaxies: starburst 
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1 INTRODUCTION 

Statistical analyses of large samples of galaxies belonging to rich 
galaxy clusters reveal a gradual trend of the quenching of star for- 
mation in galaxies at decreasing radial distance from the centres 
of clusters. This is generally interpreted as the overwhelming in- 
fluence of the gravitational field of the dark matter associated with 
the cluster, added to the effect of the hot gas in the intracluster 
medium (ICM), that serves to progressively deprive the galaxies of 
the fuel for star formation through tidal and other stripping pro- 
cesses. Thus, the core of a cluster is generally thought to be the site 
of the most dramatic transformation in the evolution of a galaxy. As 
a result, the central regions of galaxy clusters end up being domi- 
nated by passively evolving red elliptical/lenticular galaxies, while 
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their blue star-forming counterparts prefer the sparsely occupied 
regions away from clusters. 

However, over the last few years, a growing number of galax- 
ies, with an unusually high rate of star formation, have been dis- 
covered on the outskirts of rich galaxy clusters. These galax- 
ies, often spotted close to or outside the virial radius of clus- 
ters, are seen to possess tidal features, and signs of morpholog- 
ical distortions that appear to be induced by close interactions. 
A majority of thes e galaxies have been discovered beyond a red- 
shift of z = 0.2 jKeel|[2003:lMoranet al.||2005l : ISato & Martin] 
l2004lMarcillac et al]2007l ; IOemler et al.l2009T) . since at these red- 
shifts, the regions beyond the virial boundaries of clusters of- 
ten fit on to current wide-format imagers. Good examples of 
these are the infalling galaxy C153 in Abell 2125 (z = 0.246; 
IWang. Owen. & Ledl ow 2004), which exhibits comet-like tails of 
X-ray emitting gas, or the flagra ntly star-forming galaxies on the 
outskirts of Abell 2744 (z = 0.4: lBraglia et alj|2007l) . which seem 
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to lie on filaments feeding a set o f merging clusters, and in the 
lensing cluster CL0024+16 (z = 0.4; lMoran et alj .2005). where en- 
hanced star formation is detected in a string of early-type dwarfs 
outside the virial radius. 

It is not yet known whether such galaxies are common at lower 
redshifts, since surveys away from the cores of nearby clusters are 
not common. Infalling galaxies like NGC 4472, at 1.2 Mpc from the 
core of the Virgo cluster (z = 0.003), show evidence of extended 
streams of h ot X-ray emittin g gas, and a shock front in the X-ray 
observation dKraft et alj|201 l|). This is accompanied by star forma- 
tion in tidal str eams of stripped cold gas fro m its dwarf compan- 
ion UGC 7636 ( iLee. Richer. & McCalfeoOOh. Similar systems are 
seen elsewh ere in Virgo dAbramson et al.l201lh . and in Abell 3627 
(z = 0.016: ISun et alj|2010h . where the abundance of Ha and HII 
emission indicate rapid star formation. Some of these galaxies are 
plunging towards the cluster core on their own, while others like 
those comprising J0454 (z = 0.26) ar e merging with the cluster 
as a galaxy group dSchimeret"aIB"2 010). In the absence of system- 
atic studies, it is not clear how common such starburst galaxies are 
away from the cores of clusters, and indeed how important this phe- 
nomenon is in the grand scheme of galaxy evolution. One would 
like to know, for instance, whether their occurrence is related to the 
properties of the clusters, the structure of the surrounding cosmic 
web or the properties of the galaxies themselves. 

Earlier studies have suggested that evolution of galaxies 
in clusters is mainly driven by processes such as ram-pressure 
stripping, where a galaxy falling into cluster loses its gas con- 
tent through interactions with the hot ICM, thus being de- 
prived of star formation by the time it reaches the cluster core. 
Major mer gers between galaxies a re uncommon even in clus- 
ter cores jPe Lucia & Blaizotl 120071) . and galaxy-galaxy interac- 
tions dMooreet alj|l996h are in general thought to be unim- 
portant in cosmological models (e.g. iMaver et al.l |2006). How- 
ever, the observed ra dial gradient in star formation and other 
galaxy properties (e.g. Boselli & Gavazzi 2006; Porter et al JI2008I : 
iMahaian. Haines. & Ravchaudhurvl l2010h . and the existence of 
transitional (i.e. post-starburst) galaxy populations to much farther 
distances from the centres of clusters than their region of influ- 
ence (e.g. iGavazzi et al.ll201(il : IMahaian. Haines. & Ravchaudhurvl 
l201lh . indicates that the influence of the cluster does not tell the 
whole story. 

In particular, repeated high velocit y encounters be tween 
galaxies, often termed as "harassment' ' dMooreetalJ [l996). can 
play a special role in the evolut ion of galaxies in low and in- 
termed i ate density environment s ([Lewis et al. 2002; Haine s et al 



2006al; iPorter & Ravchaudhurvl 120071; iMossI 120061 ; iPorter et al.. 
20081 ; lOemler et alj|2009l ; ISmith et al]|20ld. among others). The 



frequency and impact of such galaxy-galaxy interaction pro- 
cesses depend upon the local galaxy density. So, in prin- 
ciple, any location on the cosm ic web, with gala xy den- 
sity above a 'critical' value (see Gom ez et al.l 120031 for in- 
stance), should encourage interactions, resulting in bursts of 
star formation, among galaxies . It has been shown else- 
where that filame n ts of galaxies jPorter & Ravchaudhurvl 1 20071; 
BoueetalJ 120081; iFadda. Biviano & Durretl 120081; IPorter et al.l 



20081 ; lEdwards. Fadda, & Fraveil l2010j; iBiviano et alj |201lh 



and outskirts of clusters dRines et alj 120051 ; ISmith et al.l I2OI0I 
Coppin et al. 201 1 ; Verdug oet alj|2012f) are such favourable sites. 

While simulations suggest that stripped HI and X-ray emit- 
ting gas should be ubi quitous and detectable in the central regions 
of mo st clusters (e.g. T onnesen & Brvanl201ol ; |\^i magami & Fuiital 
l201ll) . it is difficult to quantify systematic properties of galaxies 



residing on inter-cluster filaments on the cosmic web due to the 
small numbers involved on a single filament. However, by stacking 
a large number of filaments, in the Two-degree Field Galaxy Red- 
shift Survey (2dFGRS), it has been shown that the star formation in 
galaxies is enhanced 3-4 Mpc from the cores of the clusters, sug- 
gesting that occurrence of burst of star form ation in infalling galax- 
ies should be a comm on phenomenon jPorter & Ravchaudhurvl 
l2007l ; |Porteretalj2008h . 

Since outside the virial boundaries of clusters, the density of 
the ICM is too low to influence the incoming galaxies, it was sug- 
gested in these studies that the enhancement in star formation is 
due to high speed encounters occurring amongst infalling galaxies 
on the outskirts of galaxy clusters. Multiple starburst galaxies have 
been seen in similar single systems associated with filamen ts, such 
as in Abell 1763 (2 = 0.23; lFadda. Biviano & Durretl2008l) . Multi- 
wavelength studies of photometric properties of galaxies belonging 
to the Shapley supercluster (z = 0.048) revealed that the super- 
cluster environment might be solely responsible for the evolution 
of faint (> M* + 2) galaxies faaines et alfeoill) . 

This paper addresses the issue of enhancement in star forma- 
tion activity in galaxies found in the infall regions of nearby galaxy 
clusters. The star formation activity is quantified using the inte- 
grated SFR and specific star formation rate (SFR/A/*) of galaxies 
in a representative sample of local clusters, and emission line ratios 
to separate out the emission dominated by the active galactic nu- 
clei (AGN) from that of star formation. A catalogue of large-scale 
inter-cluster filaments of galaxies, sourced from the SDSS, is used 
to evaluate the impact of the position of a galaxy in the cosmic web 
on its star formation properties. 

The outline of the paper is as follows: the data and sample 
properties are elaborated in the following section. The star for- 
mation properties of cluster galaxies are described in Sj3] briefly 
discussing the behaviour of galaxies falling in the cluster for the 
first time. The impact of the immediate and the large-scale struc- 
ture (LSS) environment on the star formation properties of galax- 
ies is discussed in ij4] and the main results are summarised in Sj5] 
Throughout this work a cosmology with Q m — 1 & Qa = and 
h = 70 km s _1 Mpc -1 was adopted for calculating distances and 
absolute magnitudes. We note that in the redshift range chosen for 
this work (0.02 ^ z ^ 0.15), the results are insensitive to the 
choice of cosmology. 



2 THE SAMPLE 

We us ed the spectroscopic galaxy d ata provided by the SDSS 
DR4 dAdelman-McCarthv et all [2006), which uses two fibre-fed 
double spectrographs, covering a wavelength range of 3800-9200 
A. The resolution A/ A A varies between 1850 and 2200 in differ- 
ent bands. The galaxy magnitudes and the corresponding galactic 
extinction values in the r-band were taken from the New York Uni- 
versity Value added galaxy Catalogue (NYU-VAGC; B lan t on et al 



2005a). The fc-corrections (to z = 0.1) provided by lYang et al 



( 2007) were used to correct for the effect of distance on the mag- 
nitude of galaxies. Only those galaxies with absolute magnitude 
M r —20.5 (SDSS magnitude limit at the median redshift of our 
sample z = 0.1) were considered. 



2.1 Cluster membership 

All Abell clusters dAbell et al.lll989h in the redshift range 0.02 < 
z ^ 0.15 were considered. In order to assign galaxies to clusters 
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Figure 1. The distribution of the errors on the SFRs of galaxies in our sam- 
ple. The solid points denote the mean 68% confidence interval (C.I.) of 
the likelihood distributions for the SFRs in bins of log SFR (left) and log 
SFR/Af* (right). The error bars show the 2cr scatter in each bin. Galaxies 
with high SFR have higher S/N ratio in the emission lines, and hence tighter 
constraints on the estimated SFR and SFR/M* relative to the galaxies with 
low or no optical emission. It is also noticeable that the scatter in the C.I. in 
all the bins of SFR (and log SFR/Af *) is comparable. 



differential redshift slices were used, i.e., galaxies, within a pro- 
jected radius, R — 6 Mpc of the cluster centre at the redshift of 
the cluster, were included if they were found to be within a redshift 
slice of: 



3000 km s~ 
2100 km s" 
1500 km s" 



if R < 1 Mpc 

if 1 < R < 3 Mpc 

if 3 < R < 6 Mpc 



(1) 



This differential velocity slice approach helps in reducing in- 
terlopers without following a rigorous procedure for assigning clus- 
ter membership, and hence is well suited for a statistical study like 
this. Note that popular methods, such as the 3a v velocity disper- 
sion cut, which is often employed to define galaxy membership of 
clusters, is not our preferred choice because of our need to include 
the outskirts of clusters. 



2.2 Cluster parameters 

Once the galaxy members of the clusters were found as described 
in i j2.ll their velocity disp ersion (a v ) is determined using the bi- 
weight statistics (ROSTAT. iBeers. Flvnn & Gebhardt| [T990). Then 
assuming cluster mass M (r) oc r, r2oo (the radius at which the 
mean interior over-density in a sphere of radius r is 200 times the 
critical density of the Universe) wa s calcu lated, using the relation 
given bv lCarlberg, Yee & Ellingsor] dl997t) . 



f200 



V3 a 



(2) 



10 H(z) ■ 

Any comparison clusters (described below) which were not fully 
covered in DR 4 were excluded along with any clusters for which 
reliable velocity dispersion information was neither available nor 
obtainable from the given data. 

The final sample consists of 107 Abell clusters dAbell et all 

Il989f) with z < 0.15 lying in those regions of the sky where SDSS 
DR 4 has more than 70% spectroscopic coverage. 

2.3 Galaxy data 

The galaxy data (SFR, specific star formation ra te, SFR/M*) used 
in this paper were sourced from the catalogue of Brin chmann et al.l 



(2004, B04 henceforth) which is based on SDSS DR4. The 
Bayesian technique employed by B04 for modelling various phys- 
ical parameters including SFR and SFR/ M* is advantageous over 
using a single emission line such as Ha or [Oil] as a star formation 
indicator because the integrated SFR is a synthesized representa- 
tion of various line indices. This reduces the uncertainty inherent in 
using a single index by several orders of magnitude. The simultane- 
ous use of different indices enables us to separate the galaxies with 
a significant AGN component from those with emission entirely 
due to star formation (see below). Using the entire spectrum also 
makes it possible to obtain an independent estimate of the stellar 
mass M* (from the z-band photometry), which further constraints 
the estimated SFR/ M*, making it a more reliable parameter for 
evaluating star formation properties of galaxies, relati ve to compa- 
rable luminosity based parameter used elsewhere (e.g. iLewis et al] 
l2002t) . 

B04 divided all galaxies into star forming, AGN a nd com- 
posite using the BPT (Baldwin, Phillips & Terlev ichll98lh diagram 
and then derived the SFR and SFR/Af* from the spectra taken by 
the (3" diameter) SDSS fibre. These SFRs and SFR/ A/* were then 
corrected for aperture biases using galaxy colours measured 'glob- 
ally', and within the fibre (see B04 for details). Out of the three 
statistical estimates (mean, median and mode) for the probability 
distribution of SFR and SFR/ M* derived for each galaxy, we used 
the median of the probability distribution function (PDF) of the star 
formation parameters since it is independent of binning. 

Fig. Q]shows the 68% confidence intervals (C.I.) for log SFR 
estimate for each galaxy in our sample and the scatter in C.I. in bins 
of log SFR and log SFR/ A/* respectively. The trends are similar to 
those seen in fig. 14 of B04 (for all galaxies in SDSS DR 4). Since 
the SFR for the star-forming galaxies are measured directly from 
the high S/N emission lines (> 3 in all 4 lines: Ha, H/3, [OIII] and 
[Nil] required for the BPT diagram), the constraint on the derived 
value of SFR is tight. However, for the low S/N and non-emission 
line galaxies, most of which are passively evolving, SFR is mea- 
sured indirectly from the 4000 A break-SFR relation obtained for 
the high S/N galaxies (see B04 for details), resulting in high uncer- 
tainties in the estimated SFR and SFR/M*. We chose not to select 
galaxies based on the confidence interval values derived from the 
probability distribution function of the SFR and SFR/ Al* parame- 
ters of individual galaxies because excluding galaxies with higher 
uncertainties in SFR/M* (or SFR) distribution would exclude a sig- 
nificant number of passive galaxies from the sample. But precau- 
tionary measures are taken into account where necessary ( jj2jj) , so 
that our results remain virtually unaffected by the uncertainties in 
the SFR/ M* (or SFR) estimates for the unclassified galaxies. 

2.4 Unclassified galaxies 

Any galaxy which can not be classified using the BPT diagram, 
because of low S/N ratio or non-detection of the requisite emis- 
sion lines, is termed 'unclassified' in B04, which makes this class 
of galaxies rather heterogeneous. The spectra obtained from the 
3" fibre, which encompasses ~ 40-60% of the galaxy's light (at 
z ~ 0.1), may hence only represent the passive core of the galaxy, 
lacking any emission lines. A careful analysis of the images of 
these galaxies revealed that these are predominantly bulge domi- 
nated systems, in fact even more concentrated than AGN or com- 
posite galaxies (J. Brinchmann, private communication). A detailed 
investigation of whether all of these galaxies are indeed passively 
evolving or comprise a sub-population of dust obscured galaxies 
with an AGN supported star formation in the galactic nucleus is 
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Figure 2. (top) SFR and (bottom) SFR/Af* as a function of r-band magnitude for all the galaxies in SDSS DR4, divided according to their position in the 
BPT diagram (see £12. 31 . The contours for star-forming+Composite class (right panel) are repeated in other panels to compare the relative distributions. The 
outermost contour in each case represents 100 galaxies per bin, increasing by 100 for every consecutive contour inwards. As expected, the AGN hosts follow 
the passive galaxy distribution. It is also interesting that the difference between various classes become more distinct for SFR/M* (bottom panel). These plots 
show that absolute SFR follows the galaxy luminosity, but the trend reverses when SFR/Af* is considered. Another interesting feature that emerges here is 
that the observed bimodality in SFR seems to largely consist of the 'non-emission line' galaxies. This plot gives no indication of enhanced (or suppressed) star 
formation in AGN hosts, although this might be a result of the obscuration of emission lines by dust around the active nuclei. 



beyon d the scope of this work (but see lMahaian & Ravchaudhurvl 
2009, where we briefly address related issues). 

Fig.[2]shows the distribution of absolute SFR and log SFR/M* 
of all galaxies (-23 «S M r < -18) in SDSS DR4 as a function 
of their luminosity for the three different classes: star forming (& 
composite), AGN and the unclassified (B04). As previously specu- 
lated, the unclassified galaxy class seems to be dominated by pas- 
sive galaxies. It is also not surprising that the galaxies hosting an 
AGN which are known to be bulge-dominated systems, follow the 
same distribution as that of the unclassified galaxies. 

In agreement with previous findings, Fig. [2] shows that 
most low-luminosity galaxies are efficiently forming stars 



(Kauffman n et alj|2004 lHaines et alj|2006al. 120071 ; iDamen et all 
2008 ; lMahaian. Haines. & Ravchaudhurj201(]ir Interestingly, the 
observed bi-modality in the distribution of SFR is exclusive to the 
unclassified galaxies. If the SFR of the star-forming galaxies is a 
function of environment, one should expect the distribution of their 
SFR to show a bimodal distribution in a dataset which samples a 
wide range of environments, such as that employed here. However, 
no conclusive remarks can be made on the issue of starburst-AGN 
connection based on the analysis presented here, because at visi- 
ble wavelengths, the circum-nuclear star formation (which largely 
contributes to the 3" SDSS fibre spectra) may be entirely obscured 
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Figure 3. Top: Distribution of all (black points) cluster galaxies and star- 
bursts (red stars) within 3r200 of the cluster centres in the SFR-Af r space 
shows that the starburst galaxies span the entire magnitude range as spanned 
by other cluster galaxies. Bottom: The redshift distribution of galaxies in the 
CN (solid grey) and SB (open red) cluster samples. The K-S test does not 
distinguish between the parent population of these two samples. The red 
hatched distribution is for the starburst galaxies in the SB sample. 



by dust enshroudin g the active nucleus (e.g. IPopescu et al.ll20o3 ; 
iPrescott et al1l2007l) . 



2.5 Starburst galaxies 

A galaxy having current SFR much higher than that averaged over 
its past lifetime may be regarded as a starburst. In the context of this 
paper, dual criteria were adopted to define a starburst galaxy: log 
SFR/M* > -10 yr _1 and SFR > 10 Meyr" 1 . These thresholds 
are purely empirical, driven by the idea of selecting galaxies having 
current SFR much higher than that averaged over their past star 
formation history. These limits thus select galaxies contributing to 
the high end tail of the SFR and SFR/M* distributions only, as 
shown in Figure[3](top panel). 

In particular, they label ~ 3% of all and ~ 13% of star- 
forming galaxies in our sample of cluster gal axies as a starburst . 
In agreement with several other studies (e.g. iFeulner et al]|2005l : 
lNoeskeetai]|2007r) . we note that by choosing starburst galaxies 
according to absolute SFR or SFR/M* alone could bias the se- 
lected galaxies to giant or dwarf galaxies respectively (because 
dwarf galaxies have higher SFR/M * than giants but lower abso- 
lute SFR; e.g. iDamen et al.ll2008l) . These starburst galaxies also 
cover a wide range in structural properties. 

We stress that the above definition of starburst galaxies is used 
only to sub-classify galaxy clusters in our sample. The results pre- 
sented in this paper remain unaffected by the specific choice of 
SFR and SFR/M* thresholds quantitatively (see !j3]for a detailed 
discussion). 
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Figure 4. The distribution of red and blue galaxies (colour-selected in the 
colour-magnitude space; see text) in and around the rich cluster Abell 1 190, 
which is part of our SB sample. The starburst galaxies in its vicinity are 
shown as cyan stars and the corresponding optical images are shown inset. 
The dotted circles represent regions encompassing 1 and 3 Mpc from the 
centre of the cluster (green cross). 




D n 4000 

Figure 5. All galaxies in the SB clusters are plotted as black dots in the 
D n 4000-H<5 space. Overplotted as red stars are starburst galaxies found at 
1 ^rfeoo ^ 2. The blue circles are galaxies at the same distance from 
the cluster, but having 1 ^ SFR < 10 M0 yr -1 . Most of these infalling 
galaxies with high SFR (and SFR/M*) occupy the area in the D n 4000- 
H<5 space where galaxi es with curren tly ongoing or recent starbursts are 
expected to lie jKauffmann et al. 2003). 



3 STAR FORMATION AND ENVIRONMENT 

From a sample of 107 galaxy clusters, those having at least one 
(Mr ^ —20.5) starburst galaxy (as defined in ^2.5t . anywhere 
within 3 Mpc of the cluster centre, were selected, yielding 46 clus- 
ters termed as 'SB clusters'. The distance range of 3 Mpc was cho- 
sen since, for most of the galaxy clusters in our sample, r2oo ~ 1-5 
Mpc. Hence the starburst galaxies chosen to define this sub-class 
of clusters would either belong to the cluster or to its outskirts. In 
order to make the selection procedure independent of the uncertain- 
ties in the SFR estimate ( $2.3\ Fig.QJ, starburst galaxies which are 
unclassified ( i]2.4t but have estimated SFR ^ 10 Mq yr -1 in the 
B04 catalogue, were ignored. Redshift distribution of cluster galax- 
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Table 1. Abell clusters with starburst 
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a IPopesso eta?] f2007) 
b IPopesso et al] l2007t) 

ies in the two samples are shown in Figure [3] (bottom panel). The 
Kolmogorov-Smirnov (KS) statistic suggests that the two samples 
are statistically similar. The distribution of redshifts for the star- 
burst galaxies is also statistically similar to that of other cluster 
galaxies. 

A typical SB cluster, Abell 1190, is shown in Fig. [4] as an 
example. In this figure, galaxies are labelled as red or blue us- 
ing a colour-magnitude relation fitted to galaxies within 3Mpc. Of 
the rest of the 61 clusters which do not have any starburst galaxy 



Table 2. Comparison (CN) clusters 



Abell RA Dec Richness z Velocity 

(J2000) (J2000) dispersion 

(a v ; km s _1 ) 
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-9 51 30 
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354 
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8 1 4.0 
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0.0743 


690 


1038 


10 32 59.4 


+ 2 15 17 


1 


0.1246 


297 


1064 


10 38 46.9 


+ 1 164 


1 


0.1320 


594 


1068 


10 40 50.2 


+39 57 1 


1 


0.1386 


1317 


1080 


10 43 58.8 


+ 150 





0.1205 


432 


1139 


10 58 4.9 


+ 1 30 45 





0.0395 


249 


1169 


11 8 10.0 


+43 56 34 


1 


0.0587 


591 


1171 


11 7 29.2 


+ 2 56 31 





0.0757 


159^, 


1173 


11 9 14.4 


+41 34 33 


1 


0.0759 


570* 


1205 


11 13 23.0 


+ 2 30 29 


1 


0.0759 


810 


1225 


11 21 24.5 


+53 46 22 





0.1040 


732 


1270 


11 29 33.3 


+54 4 20 





0.0686 


600 


1302 


11 33 30.4 


+66 25 18 


2 


0.1160 


834 


1318 


11 36 30.6 


+54 58 16 


1 


0.0564 


366 


1361 


11 43 48.3 


+46 21 17 


1 


0.1167 


453 


1364 


11 43 39.8 


- 1 46 36 


1 


0.1063 


558 


1368 


11 45 1.7 


+51 15 13 


1 


0.1291 


762 


1372 


1 1 45 29.3 


+11 31 13 


1 


0.1126 


435 


1383 


11 48 13.5 


+54 37 12 


1 


0.0600 


426 


1385 


11 48 5.5 


+11 33 16 


1 


0.0831 


531 


1387 


11 48 54.1 


+51 37 12 


1 


0.1320 


753 


1390 


11 49 35.1 


+12 15 15 





0.0829 


342 


1402 


11 52 37.4 


+60 25 15 





0.1058 


219 


1452 


12 3 42.6 


+51 44 14 





0.0628 


513^ 


1552 


12 29 50.6 


+11 44 29 


1 


0.0859 


71~ 


1559 


12 33 10.0 


+67 7 31 


1 


0.1066 


693 


1564 


12 34 57.3 


+ 1 51 32 





0.0792 


687 


1566 


12 35 5.1 


+64 23 32 


2 


0.1000 


561 


1577 


12 37 52.0 


-0 16 22 


1 


0.1405 


366 


1621 


12 48 45.9 


+62 41 46 


1 


0.1029 


534 


1630 


12 5144.8 


+ 4 34 49 


1 


0.0648 


426 


1637 


12 53 59.3 


+50 49 51 


1 


0.1270 


564 


1646 


12 55 48.6 


+62 9 53 





0.1063 


468 


1650 


12 58 46.1 


- 1 45 57 


2 


0.0845 


513 


1663 


13 2 46.4 


-2 3149 


1 


0.0847 


891 


1764 


13 34 43.2 


+59 55 51 





0.1167 


585 


1773 


13 42 8.5 


+ 2 15 7 


1 


0.0766 


1098 


1780 


13 44 38.1 


+ 2 53 12 


1 


0.0786 


426 


1882 


14 14 39.8 


- 20 33 


3 


0.1367 


612 


1885 


14 13 47.6 


+43 40 15 


1 


0.0890 


1038 


1918 


14 25 9.2 


+63 9 41 


3 


0.1394 


1008 


1920 


14 27 17.5 


+55 46 46 


2 


0.1310 


570 


1936 


14 34 28.8 


+54 50 5 


1 


0.1386 


498 


2018 


15 1 12.8 


+47 16 26 


1 


0.0872 


1029 


2046 


15 12 41.4 


+34 51 9 


1 


0.1489 


948 


2062 


15 21 19.8 


+32 5 38 


1 


0.1155 


537 


2110 


15 39 43.5 


+30 42 45 


1 


0.0980 


654 


2175 


16 20 23.1 


+29 55 18 


1 


0.0972 


876 


2196 


16 27 21.6 


+41 29 43 





0.1332 


687 


2211 


16 34 4.0 


+40 57 9 


1 


0.1355 


537 


2245 


17 2 45.0 


+33 33 13 


1 


0.0852 


984 



Enhanced star formation on the outskirts of galaxy clusters 7 



within 3 Mpc of the cluster centre form the 'comparison' (CN) clus- 
ter sample. The fundamental properties of our cluster samples are 
listed in Tables [7] and [2] respectively. In total 3,966 galaxies are 
found within 3r2oo of the cluster centres, of which 3,903 are non- 
AGN and were used for all further analysis, unless stated otherwise. 
If the occurrence of starburst galaxies in some clusters is purely co- 
incidental, or the result of line-of-sight superposition, then CN and 
SB cluster samples should display similar properties. As shown be- 
low, this hypothesis fails to hold. 



3.1 The last 'bursts' of star formation 

The 4000 A break (D„4000) is a prominent discontinuity occur- 
ring in the continuum of the visible spectrum of a galaxy. It arises 
because a large number of absorption lines, mainly from ionised 
metals accumulate in a small wavelength range. In hot young stars, 
the metals are multiply ionised and hence the opacity is low, result- 
ing in a smaller D n 4000 relative to old, metal rich stars. We used 
the same definition of D n 4000 as adopted by B04, which is the 
ratio of average flux density in the narrow continuum bands (3850- 

3950 and 4000-4100 A). This definition was initially introduced by 

I 1 J 7\~. 

Balog h et alj (1199 9) and is an improvement over the original def- 
inition of iBruzuall dl983f) . in the sense that the narrow continuum 
bands make the index considerably less sensitive to reddening ef- 
fects. 

The H5 absorption line occurs when the visible light of a 
galaxy mainly comes from late-B and early-F stars. A large ab- 
sorption in YL8 implies that a galaxy experienced a major burst of 
star formation ~ 0.1-1 Gyr ago, while a stronger D n 4000 is an in- 
dicator of metal-rich ISM. Hence together EW(H<5) and D n 4000 
indicate whether a galaxy has been forming stars continuously 
(intermediate/high D„4000 and low EW(H<5)) or in bursts (low 
D n 40 00 and high EW(H<5)) over the past Gyr dKauffmann et al.l 
l2003l) . Fig.[5]shows the distribution of spectacular starburst galax- 
ies (SFR^ 10 M yr" 1 & log SFR/A/* > -10.5 yr" 1 ), rela- 
tive to all other galaxies in the SB cluster sample, in the D„4000- 
H<5 plane. It is clear that our adopted criteria for defining starburst 
galaxies on the basis of integrated SFR and SFR/A/* segregates 
well these galaxies in the space defined by spectral indices sensitive 
to recent star formation. For a comparison, we also show galaxies 
which have the same SFR/A/* but absolute SFR which is upto 10 
times lower than 'starburst' galaxies. For the latter, the distribution 
shows a larger scatter. 



3.2 Starburst galaxies on the outskirts of clusters 

We began this section by sub-classifying clusters on the basis of the 
presence of starburst galaxies in them. In order to test the hypothe- 
sis that these starburst galaxies might occur randomly within a clus- 
ter, we plot the first, square root of the second and the third moment 
(mean, standard deviation & skewness) of SFR respectively, as a 
function of (scaled) cluster-centric radius, in Fig. [6] We extend our 
analysis upto 3 r2oo bec ause it encompasses the core of the galaxy 
cluster and its outskirts jDiaferio & Gelled 19971 ; [Lewis et al.l2 002 ; 
iGomez et alj|2003l) . providing a range of galactic environments to 
explore. 

The SFR profiles of the two samples in Fig. [6] are strikingly 
different, with SB clusters showing a higher mean SFR relative 
to CN clusters at almost all radius. The striking feature of this 
stacked SFR profile is the sudden increase in the standard devia- 
tion in SFR, which, accompanied by positive skewness, indicates 



an enhancement in the mean star formation of galaxies in the SB 
clusters. This shows that several highly star-forming galaxies are 
present just outside the boundary (~1.5r2oo) of SB clusters. Given 
that the selection criteria demanded that starburst galaxies could oc- 
cur in these clusters anywhere within the given threshold (3 Mpc), 
it seems highly unlikely that occurrence of such a 'blip' is a ran- 
dom chance. There seems to be another region of increase in the 
mean SFR of galaxies in the SB clusters at ~ 0.5 r2oo- We note 
that the latter feature, closer to the cluster core, is neither due to the 
presence of group-like structures (a v ^ 300 km s _1 ; tableQJ, nor 
to the nearby interacting cluster pair Abell 2197 and 2199 present 
in our sample. Both these features in the SFR profile of SB clusters 
are thus considered real. We discuss them in detail in the following 
sections. 



3.3 Clusters & Filaments 

In order to explore whether the enhancement in SFR of galax- 
ies found on cluster outskirts ( §3.2\ is related to the properties 
of the clusters themselves, the distribution of velocity dispersion 
(<t„) and X-ray luminosities (Lx), where available, for the two 
cluster samples are plotted in Fig. |7j Since a v are derived from 
spectroscopic redshifts using an iterative procedure (ROSTAT; 
iBeers. Flvnn & Gebhardti ri990). the error in individual a v values 
is very small. Lx are obtained for 31/46 SB and 45/61 CN clus- 
ters respectively from BAX (http://bax.ast.obs-mip.fr/!, which is an 
online database providing observational X-ray measurements for 
clusters of galaxies observed by various space and ground-based 
X-ray missions. In the following we use Lx measured in the 1.5- 

2.4 keV energy band as a proxy for cluster mass, and a v as a proxy 
for the dynamical state of the cluster. 

Fig. UJ shows that despite the fact that cluster samples consist 
of a wide range of objects (group-like structures to massive clus- 
ters), their a v distributions are significantly different. As suggested 
by the KS statistic, SB clusters have higher a v than CN clusters, but 
there are no statistical differences in the distributions of the Lx for 
the two samples. To confirm that we are not selecting only the X-ray 
bright clusters in either of the samples, we further checked that the 
a v distributions of clusters from the two samples, for which Lx are 
not available, are not statistical different, as confirmed by the KS 
test. This implies that there is no difference in the mass distribution 
of the two samples. However, similar Lx but higher a v suggests 
that the SB clusters are likely to be dynamically unrelaxed relative 
to the CN sample. The observed difference in the a v distribution is 
not due to the two sample of clusters being sampled from different 
ranges of redshift since both the cluster samples cover a redshift 
range (0.02-0.15) spanning only ~ 1.5 Gyr in lookback time. 

If a cluster is being assembled via galaxies falling in through 
large-scale filaments, the unrelaxed dynamical state of the clus- 
ter will be observable as high a v , especially if one or more 
such filaments align along the line of sight. Hence the above 
observations could be interpreted as the manifestation of sys- 
tematic differences in the dynamical state of the two cluster 
samples with similar mass distributions. Such an impact of the 
global environment, in which a galaxy is embedded, on its 
star formation propertie s , is in e vitable (e.g. Haines et alj 



2011 



2008 



iPorter & Ravchaudhurvl 120071; iFadda. Biviano & Purred 
iPorter et al . 2008; Ed wards. Fadda. & FraverluOlOal) . However, an 
interesting question is whether the large-scale structure directly in- 
fluences the evolution of galaxies therein, or has an indirect effect. 
Here we use the position of galaxies along the cosmic-web to ex- 
plore this further. The catalogue of galaxy filaments used in this 
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Figure 6. The mean (thicft line, bottom panel), skewness (upper panel) and standard deviation (thin red, bottom panel) in SFR of galaxies in SB (7e/f) and 
CN (right) clusters in running bins (250 galaxies per bin) of (scaled) cluster-centric radius. Unclassified galaxies with SFR> 10 Mq yr and AGN were 
excluded. As expected, the SFR falls toward the cluster centre. But the difference in the SFR profiles of the two samples is striking. The high standard deviation 
at the cluster boundary (r ~ r2oo) indicates the presence of galaxies with extreme SFRs. It is also noticeable that the mean SFR of CN clusters is lower than 
that of SB clusters in general. 



pa per is sourced from SPSS DR4 by using th e algorithm described Table 3. Type of filaments 
bv lPimbblet. Drinkwater & Hawkrigd (2004, PDH04 thereafter). 



Just as in PDH04, all filaments are classified into morpho- 
logical classes from types I to V i.e. straight, warped, sheet-like, 
uniform and irregular distribution of galaxies respectively, with re- 
spect to an inter-cluster axis. The algorithm is well described in 
PDH04, but for the benefit of reader we briefly describe it here: to 
identify galaxy filaments in a large-scale survey such as the SDSS 
involved setting up orthogonal planar pair projections of the fila- 
ments viewed along the inter-cluster axis. The planes corresponded 
to ±20 Mpc from the inter-cluster axis and had a depth of ±5 Mpc 
around it. This choice of depth and width ensures that the ability 
to visualise highly curved and complex filaments is retained (fig. 1 
of PDH04 illustrated this process in more detail). Finally, those fil- 
aments for which both Mahajan and Pimbblet agreed on the same 
morphological class were included in the 'clean' sample, and were 
used in all further analysis. 

Simulations have shown that properties, such as matter density 
of filaments, are closely related to their morphology dColberg et al.l 
l2005l ; lDolag et alj2006h . Therefore, by studying the morphology of 
filaments through which a galaxy approaches the cluster, we intend 
to explore how, if at all, is the LSS modulating galaxy evolution. 
Fig. [8] shows the distribution of morphological type of filaments 
belonging to both the cluster samples. The errors were calculated 
using binomial statistics, where, if within a fixed radial range, n out 
of N clusters in a given sample had a particular type of filament 
feeding them, then the likelihood of a clus ter to have that class o f 
filament is given by L oc / n (l - f l ) N ~ n jp e Propris et al.l2004l) . 



Sample 


I 


II 


III 


IV 


V 


undetermined 


SB 


18 


10 


2 


1 


7 


47 


CN 


25 


35 


6 


3 


9 


68 



The error on the fraction of clusters is given by the standard devia- 
tion 



A/, = 



n(N - n) 



JV :i 



(3) 



Each distribution was normalised to unity. The number of filaments 
classified for each cluster sample are listed in Table [3] The choice 
of the magnitude cut of M r ^ —20.5 excludes all low luminosity 
galaxies from the analysis, thus significantly underestimating the 
density of some of the galaxy filaments. Hence the morphological 
class for a large fraction of the galaxy filaments could not be de- 
termined by visual inspection of density contours on the projected 
planes. On the other hand this also implies that the uncertainty in 
determining the morphological class of the 'clean' filament sample 
is very small. 

Fig. [8] suggests that the CN clusters are more likely to be fed 
by Type II (warped) filaments. The number of filaments associated 
with individual clusters do not show any particular trend. 
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Figure 7. The velocity dispersion (left) and X-ray luminosity (Lx) (right) distribution of SB (red open) and CN (grey solid) cluster samples respectively. The 
Lx are available for almost 70% of both samples. While the SB clusters are found to have systematically higher velocity dispersions than the CN clusters, 
the distribution of Lx show no significant statistical differences, suggesting that the SB clusters although not more massive than their CN counterparts are 
dynamically more unrelaxed. 



4 DISCUSSION 
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Figure 8. The distribution of the morphologi cal type of filament (according 
to IPimbblet. Drinkwater & Haw krigg 1 2004) classification) for the 'clean' 
filaments in the SB (dotted) and CN (solid) cluster samples respectively. 
The errors are calculated using Eqn[5] The SB and CN clusters show signif- 
icantly different preferences for type I and type II filaments respectively. 



Superclusters are the most massive and youngest structures 
in the Universe, within which groups and clusters of galax- 
ies evolve through interactions and mergers. In this hierarchi- 
cal model, clusters are assembled by accreting surroundin g mat- 
ter in the form o f galaxies and gala xy groups (e.g. iBekkil 
1 19991 ; Istrackll200d;l0emleretalj|2009h through large-scale fila- 
ments fceldovich. Einasto & Sh~n darin 1982; Colberg elaflbOOd ; 
Porter & Ravchaudhurvl 120071; IFadda. Biviano & Durretl 120081 : 



Porter et al.l 120081 : lEdwards. Fadda. & Fravei 



On 



course of their journey from the sparsely populated field to the 
dense cores of clusters, galaxies encounter a range of different en- 
vironments, under the influence of which they evolve, not only by 
exhausting their gas content in forming stars, but also by losing the 
gas to their surroundings. Thus, one of the most crucial links re- 
quired to fully understand galaxy evolution is to find how do the 
galaxies in cluster cores become passive while those in the field 
remain star-forming. 



4.1 Excessive star formation on the outskirts of clusters 



Studies to understand the evolution of galaxies in vari- 
ous environments generally involve large-angle sky surveys 
(e.g. [k auffmann et al. 200 4]), or target particular environments 
dGavazzi et ail 120031: iBalogh et alj |2004: iMercurio et all |2004|; 



iRines et alj 120051 : IVerdugo et al] 120081) . Despite the well-known 



dominance of passively evolving galaxies in clusters, several 
starburst galaxies have been serendipitously discover ed near the 
boundaries of intermediate redshift cluster s (z ^ 0.2; Keel 2005; 
Moran et al] 120051; ISato & Martini 120061; iMarcillac et al] 



2007 



Fadda. Biviano & Purred 120081; lOemler et al] 120091), as well as 
in the local Universe iSun & Murray! |2002| ; iGavazzi et afl 120031 : 
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Figure 9. True colour SDSS images of some of the typical starburst galax- 
ies found on the outskirts (l-2r20o) of SB clusters in g', r' and i' filters 
(represented by red, green and blue colours respectively). These galaxies 
cover a wide range in morphologies, and some show structural deformities 
and presence of close neighbours. The star-forming regions are predomi- 
nantly nuclear, but this might not be surprising since the star formation rate 
is based on fibre spectra. 

ICortese et ai]|2007l ; iReverte et al1l2007t) . An interesting question 
then is whether the occurrence of such starburst galaxies in some 
clusters is due to chance, or do they constitute a missing link in un- 
derstanding galaxy evolution and the assembly of the cosmic web. 

In a study of inter-c luster filamen ts derived f rom th e 2dFGRS, 
IPorter & Raychaudhury! d2007t> and IPorter etal] ( 120081) found that 
the mean SFR of galaxies (z < 0.1) in inter-cluster filaments in- 
creases at ~ 3 Mpc from the cluster centre. However, due to the un- 
availability of requisite data for separating AGN from star-forming 
galaxies, and of SFRs derived using sophisticated algorithm^ their 
results had some inherent uncertainties. Their results found support 
in observations of extended fields around lensing clusters, where 
iMoran et al.l J2005I) , for example, had found that all of their [Oil] 
emitting galaxies in CL0024 (z — 0.4) lie around t he virial ra- 
dius o f the cluster. Considering a simple infall model, Moran et al. 
(2005) deduced that the observed excess in the [Oil] emitters at 
the cluster's periphery can be explained in a scenario whereby an 
early-type galaxy on its infall experiences a burst of star forma- 
tion at the virial radius lasting 200 Myr and contributing 1% of 
the total galaxy mass. They attributed the burst of star formation 
to interaction of galaxies with shocks in the ICM or the onset of 
galaxy-galaxy harassment. 

Using the SDSS a nd 2dFGRS spectroscopic data for galax- 
ies, Balogh et alj fc004h found a critical local projected density of 
E5 ~ 2 Mpc -2 in the SFR-density relation characterized by the 
near complete absence of galaxies with large EW(Hq) at densities 
greater than that (see their figure 3). This critical density is sim- 
ilar to that found near the cluster boundary (~ 2 — 3h 2 Mpc -2 ; 
iRines et alj|2005h . 

By analysing a sa mple of close galaxy pairs in the SDSS, 
ISolAlonso et alj j2006l) have suggested that galaxy-galaxy inter- 

1 Porter et al. (2007, 2008) used a statistical quantity, the r;-parameter from 
the 2dFGRS database, which is derived from the prin cipal component anal- 
ysis of the spectra as a proxy I Madg wick et alj2002l) . 



actions are most efficient in inducing a burst of star formation 
when the projected distance between galaxies is ~ 100 kpc and 
their relative velocity is ~350km s" 1 . Furthermore, they found 
that these conditions are most effective in low/intermediate density 
environments. The conditions just outside the cluster boundary are 
thus well suited for galaxies to tidally influence each other. How- 
ever, given the complicated kinematics involved, galaxies do not 
have enough time to merge together, as they progressively drift to- 
wards the cluster core with an additional effective infall velocity of 
~ 220-900 km s _1 . For instance, even for a broad range of cluster 
masses 10 12 < A/ C i us < 10 15 MM. IPivato. Padilla & Lambasl feOOd) 
find that the maximum infall velocity occurs in a narrow cluster- 
centric radial range of 2 < r max < 6 Mpc. 

If the dominant process is galaxy-galaxy harassment in the 
infall regions of clusters, successive high-speed encounters be- 
tween galaxies will lead to the inflow of gas towards t he central re - 
gions of the harassed (mostly low-mass) galaxies (e.g. Fujita 1998). 
lLake, Katz. & Moord ( fl998h , for instance, showed that most of the 
fuel for star formation in a harassed galaxy may be driven to the 
central kpc of a galaxy within 1-2 Gyr in such interactions, and 
in many cases in their simulations, half of the gas mass was trans- 
ferred to the core of such a galaxy in < 200 Myr. This will cause 
rapid star formation in the core, and of course, if the galaxy is 
massive enough to have an AGN, it might be fuelled as well. In 
Fig. 8, we present a montage of some of the "starburst" galaxies 
within a radial distance of 1-2 r2oo of the centres of our SB clus- 
ters, where we have used g'r'i' images from the SDSS archive and 
produced "true colour" images. We can see that these galaxies rep- 
resent a wide variety of morphologies, though they seem to be pre- 
dominantly disk galaxies, and many have close neighbours. What 
is striking from these images is that most have concentrated nu- 
clear star formation, which is indicative of such processes in action. 
One caveat is that our galaxies were chosen on the basis of spectra 
obtained with fibres, so, barring the farthest galaxies in the sam- 
ple, many spectra have preferentially sampled the central regions 
of these galaxies. 

The density of the ICM at the periphery of clusters is small. 
Hence, any cluster related environmental process that could directly 
influence the infalling galaxies can be ruled out as a cause of en- 
hancement in the SFR of galaxies (but see iKapferer et al.l I2009L 
whose simulations suggest that under some particular conditions 
the ICM density on cluster outskirts is enough for ram-pressure 
stripping of infalling galaxies). The local conditions on the out- 
skirts of clusters, i.e. small relative velocities between pairs of 
galaxies, and higher galaxy density with respect to the field, seem 
to favour interactions among galaxies, and may lead to the ob- 
served enhancement in the S FR of galaxies in these regions. Several 
other observational studies (Lewis et al. 2002; Balogh et al .1120041; 
I Moran et alj2005tlSol Alonso et alJl2006l ; |Porter & Raychau dhury 
20071; IPorter et alj 120081) and simulations lFuiital 1 19981 ; iGnedinl 
2003b) have also concluded that galaxy-galaxy interactions are 
the most suitable explanation for the observed differences between 
galaxies in and around clusters, and those in the field. 

The 'blip' near r2oo in the mean SFR profile of SB clusters 
(Fig. [6} could be a result of these short-time interactions. As a 
galaxy approaches the core of the cluster into which it is falling, it 
begins to experience the dynamical and tidal influence of the resid- 
ual substructure in the cluster ( IGnedinl l2003bl) . We would like to 
stress that the reason such a trend has not been widely noticed is 
that it is very difficult to detect and quantify the enhancement in 
the star formation activity, occurring on such small timescales and 
confined to such narrow regions on sky jPivato, Padilla & L ambas 
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Figure 10. Fraction of star-forming galaxies in the SB (red circles) and 
CN (black squares) cluster samples respectively, as a function of scaled 
cluster-centric distance. The distributions shown by solid lines correspond 
to the selection of star-forming galaxies as log SFR/Af* —10.5 yr 
& SFR ^ 10 Moyr - 1 , while for the dot-dashed lines the criteria are log 
SFR/Af* > -10. yr -1 & SFR > 1 M Q yr-\ respectively. The fraction 
of star-forming galaxies in the CN clusters is almost always lower than the 
SB clusters, irrespective of how star-forming galaxies are selected. 



2006; Porter et al ]Hxl). Thus, it is not surprising that despite 
stacking a fairly significant sample of clusters, only a subdued sig- 
nal was detected here. 

FigurefTOlprovides some more insight into the difference in the 
SFR profile of the SB and CN clusters. A comparison of the fraction 
of star-forming galaxies as a function of the cluster-centric radius 
for the two samples shows that the fraction of starburst galaxies is 
higher in the SB clusters at almost all radius, irrespective of how 
the star-forming galaxies are selected. This evidently shows that 
the CN clusters are likely to be dynamically relaxed, virialized sys- 
tems, where the total rate of star formation is lower than the SB 
clusters. 

The other peak in the mean SFR profile (Figure [6}, observed 
nearer the cluster centre, could be a consequence of the onset of the 
cluster-related processes such as ram-pressure stripping, or it could 
be the result of an alignment of several cluster-feeding filaments 
along the line-of-sight. While the confirmation of latter requires 
data at other wavelengths and sophisticated modelling, the former 
can only be tested in simulations. We leave this issue for later con- 
sideration; 

iKauffmann et aljj2004 their Fig. 12) showed that the D„4000 
feature in star-forming galaxies is well correlated with the in- 
ternal extinction (measured in SDSS z-band). For D n 4000 ~ 1 
(most of our galaxies have D„4000> 1; see Fig. [5}, the in- 
ternal extinction can be as high as 1.5 magnitudes (in the z- 
band). Hence it is very likely that the number of starburst galax- 
ies observed at visible wavelengths may be severely underes- 
timated, making observations of this phenomenon even more 
difficult with data in one wavelength regime alon e (also see 
Gavazzi et all 120031; IMahajan & Raychaudhury 2009; IWolf etal.1 



20091 ; IMahajan. Haines. & Ra vchaud hury | |2oT 



A cluster with a deep potential well attracts relatively more 
galaxies compared to one with a shallower potential, which im- 
plies that the increased galaxy density on the outskirts of a cluster 
implies higher probability for galaxy-galaxy interactions in more 
massive clusters. In the literature, however, there seems to be an 
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Figure 11. Fraction of galaxies (1.4 < r/r 2 oo =S 2.8) in the SB (doted) 
and CN (solid) cluster samples as a function of scaled relative velocity. 
No significant peak at or above implies that these ga laxies are a mix of 
infalling and backsplash galaxies jGill. Kn ebe & Gibson 20051). 



ambiguity in the observational results. In a study of 25 clusters 
(0 < z < 0.8), Poggiant il-tlii] J2006l) found that although the 
fraction of emission-line galaxies is anti-correlated with the clus- 
ter's a v at high redshift, the intermediate mass (cr„ ~ 500-600 km 
s -1 ) z ~ clusters do not show any systematic trend. lBiviano et all 
showed that at z ^ 0.1 the fraction of emission-line galaxies 
is higher by ~ 33% in clusters with a v > 900 km s _1 as compared 
to those with a v < 600 km s ~ . On the o t her ha nd, in a study of 
17 clusters from the 2dFGRS, lLewis et all J2002I) found no differ- 
ence between the mean luminosity normalised SFR (proxy for spe- 
cific star formation rate) profile of their clusters separated around 
a v — 800 km s _1 . In this paper, using Lx as a proxy for the mass 
for almost 70% of our clusters, we do not find any difference in the 
mass distribution of the two samples, but we find the distributions 
of the a v to be significantly different. We believe this is because 
a v represents the dynamical state of a cluster, i.e. the degree of re- 
laxation, rather than the mass of the cluster, as is often used in the 
literature. 

It might be argued that the observed enhancement in the 
SFR of galaxies, on the periphery of c lusters, is due to the 
mixing of 'backsplash' galaxies (e.g. Gill, Knebe & Gibson! 



120051 ; IMahajan. Mamon. & Ravchaudhurvl 201 1[ and references 
therein) with those falling into the cluster for the first time. This 
'backsplash' population comprises galaxies which have crossed 
the cluster core at lea st once and are observed on the other 
side of its periphery dBalogh et al] |2000|; iMamon et at] |2004|; 



i Rines et alj2005l:lGill. Knebe & Gibsonll2005l) . In their simulation, 
bill, Knebe & Gibson! ( |2005» " showed that backsplash galaxies can 
be detected by their distinct centrally peaked velocity dispersion 
distribution in the radial interval 1.4-2.8r2oo, while a non-zero ve- 
locity distribution peak shows that the galaxies are falling into the 
cluster for the first time. 

This hypothesis was tested by plotting the distribution of the 
relative line-of-sight (los) velocities of galaxies in the interval 1.4- 
2.8r2oo- The distribution for SB (and CN) cluster samples, as shown 
in Fig.Qj] showed no prominent peak at \Av\/cr z = 0. This im- 
plies that the galaxies in the infall regions of the cluster samples 
studied here are a mix of infall and backsplash populations. This is 
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in agr eement with earlier studies (Rin es et alj2005l ; |Pimbblet et al.l 
2006) confirming that some emission line galaxies in and around 
clusters are a result of backsplash, but a pure backsplash model 
does not explain all their properties. Using simple m odels on sim- 
ilar data, Mahaj an. Mamon. & R avchaudhurvl J201 ll) have shown 
that even a single passage through the cluster core is enough to 
completely quench the star formation in a galaxy. 

Earlier studies showed that SB galaxies in di stant (z ~ 0.5) 
cluste rs are likely to be low luminosity dwarfs (Poggia nti et al.l 
1 19991) . As mentioned before, in the observation of enhanced 
[Oil] emission at 1.8 Mpc from the centre of CL0024 (z = 
0.4; iMoran et alj l2005t) . all of the emission-line galax ies are 
dwarfs (Mb > -20). Elsewhere, lHaines et al] feOllh found 
an excess of blue dwarf galaxies (M r > —18) around 
1.5 Mpc from the centre of the Shap ley supercluster, while 
Mahajan, Haines. & Raychaudhury! d20ld) discovered a similar in- 
cidence of excessive blue dwarfs on the outskirts of the Coma clus- 
ter (z = 0.03). In the Coma supercluster (z — 0.03), dwarf galax- 
ies exhibit a stron ger star formation-density relation than the ir mas- 
sive counterparts (Mahajan, Haines, & Ravchaudhurvll2010h . In an 
acc ompanying study of dwarf galaxies in th e Coma superclus- 
ter, [Mahajan, Haines, & Ravchaudhuryj d201ll) showed that all the 
post-starburst k+A galaxies are dwarfs, and are found in or around 
the clusters and groups in the supercluster. 

These r esults are in agreement with those presented by 
ISmith etltil J2008h . who showed that the red-sequence dwarf 
(M r > —17.5) galaxies on the outskirts of the Coma cluster have 
been truncated of star formation. Due to the flux limitations of the 
SDSS data used in this work, the low luminosity (M r ^ —20.5) 
galaxies could not be included in this analysis, but the picture por- 
trayed above will hold for them too. The only exception in case of 
dwarf galaxies might be the absence of the second burst seen near 
(~ 0.5r2oo) the cluster centre (Fig. [6j. This is because amongst 
others, e vaporation of cold gas becomes very important for dwarf 
galaxies (Fujita 2004; Bosel li & Gavazzil2006T) . Hence even a rela- 
tively softer encounter with a more massive galaxy may lead to loss 
of all gas via tidal stripping and/or evaporation. 

4.2 Star formation and the large-scale structure 

In order to evaluate their role in the evolution of galaxies, it 
is crucial to understand the properties of galaxies on superclus- 
ter filaments. High resolution simulations suggest that low over- 
density regions like galaxy filaments can host a significant frac- 
tion of gas with temperatures between 10 5 -10 7 K dCen & Ostriked 



1999| ; ICroft et alj20qil;|Phillips. Ostriker & Cenl200ll ; lDolag et al.l 



20061: lEdwardsetallg OlOb). This gas may be detected as ra- 
dio emission (e.g. lKimetai]|l989l), or through bremsstrahlung 



emission in soft X-rays (e.g. 



)6J). But not only 
is the signal too weak for the present day detectors, the for- 
mer can be employed only for very nearby cluste r pairs. Hy- 
drody namical dark matter (+ gas) simulations (e.g. iDolag et all 
2006) find that due to very low gas densities and low temper- 
atures (< 10 7 K), it is very difficult to observe these fila ments 
directly. Recent Suzaku observations dMitsuishi et alj |2012f) seem 
to have confirmed (also iKull & Bohringed 1 19991 from ROSAT) 
an excess of hot X-ray gas in the filament joining the two rich 
Abell clusters in the core of the Shapley Supercluster. These 
observations are currently very difficult, and with no advanced 
X-ray observatory in the horizon, indirect ways of understand- 
ing the role of filaments in the evolution of cosmos, by study- 
ing the effect of the filamentary environment on the galaxies are 




Mean distance to the three nearest neighbours (Mpc) 

Figure 12. The mean distance to the three nearest neighbours, within a red- 
shift slice of ±2, 000 km s -1 , for galaxies l-2r2oo from the cluster centre. 
The purple distribution shows the starburst galaxies, while the grey hatched 
histogram represents other cluster galaxies at the same radial distance from 
the cluster centre. We exclude galaxies which do not have at least three 
neighbours within this range. 



necessary (e.g. | Porter & Ravchaudhurvl 120071; iBoue et alj I2OO8I ; 
iFadda, Biviano & Durretll2008l : |Porter et alj|2008l). 

By employing semi-analytic modelling, Gonzal ez & Padillal 

(2009) showed that galaxies with similar local density but dif- 
ferent 'global' environment (i.e. different positions on the cos- 
mic web), have subtle but non-negligible differences in their 
SFR and colours. They attribute this result to the coherent mo- 
tion of galaxies approaching the cluste r s through large-scale fil- 
aments ( Porter & Ravchaudhurv 2007; Fadda, Biviano & Durret 



l2008l:IPorter et alJl2008l : lEdwards. Fadda. & Fravei2010al) . IStruck 
d2006l) discussed an interesting scenario, whereby galaxies in a 
small galaxy group crossing a cluster may be pulled towards each 
other, thus increasing the local galaxy density by an order of mag- 
nitude, and the probability of galaxy-galaxy interactions by about a 
factor of 100 (density squared). Fig. 1 121 pro vides some evidence in 
support of this sc enario (also see lGavazzietal .1 2003; Mossl l2006l: 
lOemler et ai]|2009l) . For galaxies at 1 < T2oo =S 2 from the cluster 
centre, the mean projected distance to the three nearest neighbours 
within ±2, 000 km s _1 of the starburst galaxies is < 1 Mpc in 78% 
of the cases. However, for all other cluster galaxies in the same ra- 
dial range, this is true only about 66% of the time, suggesting that 
at a given cluster-centric distance, starburst galaxies are more likely 
to be in a relatively denser environment. This is worth verifying in 
larger samples as they become available. 

Type I (straight) filaments exist between close (< 5- 
10 Mpc) cluster pairs dColberg et alj 120051 ; IDolag et alj l2006t 



Pimbblet, Drinkwater & Hawkrigg 2004 ). However, the majority 
(^ 50%) of the filaments in the cosmic web are of type II (warped). 
Some attempts have also been made to unde rstand the density 
profile of type I filaments IColberg et alj|2005 | . also see Edwards 
et al. 2010b). The i r relatively shorter length dColberg et alj|2005l : 
IDolag et alj 120061 : [Pimbblet, Drinkwater & Hawkrigd 120041) and 
higher matter density may imply that the probability of galaxy- 
galaxy interactions occurring amongst galaxies on these filaments 
is higher than on their counterparts. Fig.[8]favours this scenario. It 
shows that ~ 32% of the (morphologically classifiable) filaments 
feeding the comparison (CN) clusters are of type I (straight), as 
compare d to ^47% in the SB c lusters. This is also in sync with the 
results o f IColberg et ail (2005), who showed that in the cold dark 
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matter Universe, more massive structures are more strongly clus- 
tered than their less massive counterparts. Hence, the unrelaxed SB 
clusters (a v > 500 km s _1 ) are more likely to have a straight fil- 
ament connecting them to a nearby companion cluster. We note 
that the detection of a significant fraction of our filaments remains 
uncertain because of low galaxy (M r —20.5) densities, so the 
quantities mentioned in this section may only be considered repre- 
sentative and not precise. 

On a final note, we suggest that SB clusters are most likely dy- 
namically unrelaxed clusters, presently being assembled via galax- 
ies falling in through straight filaments. The occurrence of starburst 
galaxies on the periphery of these clusters is a direct consequence 
of the way local conditions (i.e. galaxy density and relative veloc- 
ity) are affected by the global layout of the cosmic web. 



4.3 Future implications 

Compared to galaxies at z ~ 0, distant galaxies are more gas-rich. 
Hence it would not be an exaggeration to say that such events of 
spectacular starbursts would be more commonly seen at higher red- 
shifts. Confirmation of this, however, can only result from the study 
of an unbiased lis t of high redshift clusters, sampled upto the ir in- 
fall regions (e.g. Kartaltep eet alj|2008l : IVerdugo et al.ll2008l also 
see Coppin et al. 201 1 for recent observations of star formation re- 
lated FIR emission on the outskirts of clusters at z ~ 0.25). 

Once such data are available, one can ask whether the ex- 
cess of star-forming galaxies observed in clu sters at higher red- 
shifts , the well known 'Butcher-Oemler effect' ( [Butcher & Oemlerl 
1 1984b is linked to the occurrence of starbursts in infalling galax- 
ies, which appear to be nearer the cluster centre in projection (e.g. 
lOemler et alj|2009l). For a sample of X- ray-luminous clusters at 
0.18 < z < 0.55. lEllingson et alj feOOll) showed that the galaxy 
population in the cluster cores are the same at all epochs, but higher 
redshift clusters have a steeper gradient in the fraction of blue 
galaxies when the regions outside the core (> 0.5rsoo) are con- 
sidered. Such colour-based studies, however, may have inherent 
problems due to the non-negligible fraction of galaxies that devi- 
ate from the normal correspondence be tween photometric colours 
and s pectroscopically derived SFR (see M ahaian & Ravchaudhurvl 
2009, for instance). Another interesting idea to explore, though rel- 
atively more plausible in simulations, is finding the likelihood of 
detecting a starburst galaxy on a filament feeding a cluster along the 
line of sight. Such a galaxy will appear near the centre of cluster, 
yet showing physical and intrinsic properties as that of its counter- 
parts far away from the cluster core (for instance see lOemler et al.l 
l2009h . 

It would be very interesting to study the intermediate 
mass/dwarf (M r > —20.5) galaxies, in the infall regions of the 
SB clusters, which may be harassed even more easily by their rela- 
tively massive counterparts. One needs to perform a detailed mor- 
phological analysis of the starburst galaxies, such as those studied 
in this work. Considering galaxy-galaxy harassment to be the dom- 
inant environmental mechanism affecting the star formation prop- 
erties of galaxies on cluster outskirts, one expects to see signatures 
of tidal interactions with ne ighbouring galaxie s in visible or other 
short wavelength data (e.g. ISmifh et ai]|2010l) . Unfortunately, the 
spatial resolution of the SDSS photometry, used in this work, is not 
suitable for this, which we leave for future follow-up. 



5 SUMMARY AND CONCLUSIONS 

Using the integrated SFR (and SFR/M*) derived by B04 for galax- 
ies in the SDSS DR4, the star formation properties of galaxies 
(M r < -20.5) in the nearby (0.02 < z < 0.15) Abell clusters 
were analysed. The main results are summarised below: 

• Star formation activity in galaxies on the outskirts of some 
clusters is enhanced near their periphery (r2oo) and further into the 
clusters at ~ 0.5r2oo- In agreement with several other observational 
and simulations' studies, we find that this enhancement in the SFR 
of galaxies on the clusters' periphery is most likely due to enhanced 
galaxy-galaxy interactions. 

• Starburst galaxies on cluster outskirts inhabit relatively denser 
environments, compared to other cluster galaxies at similar dis- 
tances from the cluster centre. 

• Clusters with evidence of enhanced star formation activity on 
their periphery are likely to be dynamically unrelaxed, and fed by 
straight and well-defined supercluster filaments. 

• The infalling galaxies (1.4 ^ r/r2oo ^ 2.8) in our clusters are 
a mix of galaxies that have passed through the cluster core at least 
once (the 'backsplash' population), and those which are falling into 
the cluster for the first time. 
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